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Abstract 

We consider the lowest order radiative corrections for the decay K —* 
■K°e v, usually referred as K^ z decay. This decay is the best way to extract 
the value of the Vus element of the CKM matrix. The radiative corrections be- 
come crucial if one wants a precise value of V us . The existing calculations were 
performed in the late 60's Q |^ and are in disagreement. The calculation in Q 
turns out to be ultraviolet cutoff sensitive. The necessity of precise knowledge of 
V us and the contradiction between the existing results constitute the motivation 
of our paper. 

We remove the ultraviolet cutoff dependence by using A.Sirlin's prescription; 
we set it equal to the W mass. We establish the whole character of small lepton 
mass dependence based on the renormalization group approach. In this way we 
can provide a simple explanation of Kinoshita-Lee-Nauenberg cancellation of 
singularities in the lepton mass terms in the total width and pion spectrum. We 
give an explicit evaluation of the structure-dependent photon emission based on 
ChPT in the lowest order. We estimate the accuracy of our results to be at the 
level of 1%. The corrected total width is V = r (l + S) with S = 0.02 ± 0.0002. 
Using the formfactor value /+(0) = 0.9842 ± 0.0084 calculated in (m) leads to 
\V ua \ = 0.2172 ±0.0055. 
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1 Motivation 




a) b) c) 

Figure 2: real photons 



For corrections due to virtual photons see figure |l|, for corrections due to real 
photons see figure ||. 

The -Ke3 decay is important since it is the cleanest way to measure the V us matrix 
element of the CKM matrix. If one uses the current values for V u d, V us , and V u b taken 
from the PDG then |14d| 2 + |Vu S | 2 + \V u b\ 2 misses unity by 2.2 standard deviations 
which contradicts the unitarity of the CKM matrix and might indicate physics beyond 
the Standard Model. The uncertainty brought to the above expression by V us is about 
the same as uncertainty that comes from V u a- Therefore reducing the error in the V us 
matrix element would reduce substantially the error in the whole unitarity equation. 
Reliable radiative corrections, potentially of the order of a few percent are necessary 
to extract the V us matrix element from the K e ^ decay width with high precision. 
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Calculations of the radiative corrections to the decay were performed inde- 
pendently by E.S.Ginsberg and T.Becherrawy in the late 60's J^, pL Their results for 
corrections to the decay rate, Dalitz plot, pion and positron spectra disagree, in some 
places quite sharply; for example Ginsberg's correction to the decay rate is —0.45% 
while that of Becherrawy is —2% (corresponding to corrections to the total width 
T of 0.45 and 2 respectively). We have decided to perform a new calculation since 
results of the experiments will become available soon and to explore the causes of the 
discrepancies in the previous calculations. Recently a revision of E. Ginsberg's paper, 
with numerical estimation of the radiation corrections |l4j| was published. 

Our paper is organized as follows. The introduction (Section 2) is devoted to 
the short review of kinematics of the elastic decay process (without emission of real 
photon). In Section 3 we put the results concerning the virtual and soft real photons' 
emission contribution to the differential width. In Section 4 we consider the hard 
photon emission including both the inner bremsstrahlung (IB) and the structure- 
dependent (SD) contributions and derive an expression for the differential width by 
starting with the Born width and adding the known structure functions in the leading 
logarithmical approximation (the so-called Drell-Yan picture of the process). We give 
the explicit expressions for the non-leading contributions. In Section 5, we summarize 
our results and compare them with those in the previously published papers. 

Appendix A contains the details of calculations of virtual and real soft photons 
emission. 

Appendix B contains the details of description of hard photon emission both by 
IB and SD mechanisms. Our approach to study the hard photon emission differs 
technically from the ones used in papers |IJ |) . 

Appendix C contains the explicit formulae for description of SD emission including 
the interference of IB and SD amplitudes. 

Appendix D is devoted to analysis of Dalitz-plot distribution and the properties 
of Drell-Yan conversion mentioned above. 

Appendix E contains the list of the formulae used for the numerical integration. 

Appendix F contains the details of kinematics of radiative kaon decay and,besides 
the analysis of relations of our and paper 0] technical approaches. 

In tables 1,2 and graphics (fig. (3,4,5,6)) the result of numerical estimation of Born 
values and the correction to Dalitz-plot distribution and pion and positron spectra 
are given. 



2 Introduction 

The lowest order perturbation theory (PT) matrix element of the process K + (p) — ► 
7r°(p') + e + (p e ) + v(pv) has the form 

M = ^V: a F v (t)u( Pv ) lv (l + l5 )v( Pe ) (1) 
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where F u (t) = -^(p + p') v f+(t). Dalitz plot density which takes into account the 
radiative corrections (RC) of the lowest order PT is 

lH -Ca (y,z)(l + 6(y,z))(l + \*) = ^ofr. Z) (1 + S(y,z)), (2) 



dydz \ ra\ ) dydz 



dydz, r = f — — —dydz, 



dydz J dydz 

where the momentum transfer squared between kaon and pion is: 

t = (p - p'f = M 2 K {1 + r n -z) = M 2 K R(z) . 

We accept here the following form for the strong interactions induced form factor 
/+(*)= 

f+(t) = f+(0)(l + X + ^\ , (3) 

according to PDG A + = 0.0276 ± 0.0021. From now on we'll use M 2 instead of M\. 
We define 

c = M5 f^ us? \um\ (4) 

and 

a (y, z) = {z + y - 1)(1 - y) - rv + 0(r e ) . (5) 
Here we follow the notation of : 

r e = m 2 e /M 2 , r^=m 2 jM 2 ; (6) 

where m e , m^, and Mk are the masses of electron, pion, and kaon; two convenient 
kincmatical variables are 

y = 2pp e /M 2 , z = 2pp'/M 2 . (7) 

In the kaon's rest frame, which we'll imply throughout the paper, y and z are the 
energy fractions of the positron and pion: 

y = 2E e /M, z = 2E n /M . (8) 

The region of y, z-plane where ao(y, z) > we will named as a region D. Later, when 
dealing with real photons we'll also use 

x = 2w/M . (9) 

with w-real photon energy. 

The physical region for y and z (further called D -region) is (J] 

2^/rl < y < 1 + r e - r v , 

Fi(y) - F 2 (y) < z < F^y) + F 2 (y) , 

F^y) - (2 - y)(l + r e + - y)/ [2(1 + r e - y)\ , 

F 2 (y) = Vy 2 - 4r e (l + r e - ry - y)/ [2(1 + r e - y)] ; (10) 
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or, equivalcntly, 



2V^V < z < 1 + r ^ - r e , 
F 3 (z) - F 4 {z) <y< F 3 (z) + F 4 {z) , 
F 3 (z) = (2 - z)(l + rv + r e - z)/[2(l + r w - z)] , 



F 4 (z) = y/z 2 -4r«(l + r w - r e - z)/[2(l + - z)] 



(11) 



For our aims we use the simplified form of physical region (omitting the terms of the 
order of r e ): 

2^ < V < 1 - r„ , c(y) < z < 1 + »y (12) 

with 

c{y) = l-y+ 

2V^<z<l + r v , b-{z)<y<b(z) (13) 

b-(z) = l-±(z+y/z*-4r v ) , 

Hz) = 1 - X - (z - y/z 2 - 4ry) . 
For definiteness we give here the numerical value for Born total width. It is: 

G 2 F M 5 K \V us f+(0)\ 2 f , f , * 



or, 



with 



64tt 3 



/*/ 



dza (y,z)(l + X+^rY 



rat 



5.36|K,/+(0)| 2 x lQ- 14 MeV. (14) 



Comparing this value with PDG result: (T K+e3 ) exp = (2.56 ± 0.03) x 10~ 15 MeV we 
conclude 

(K S /+(0)|) Q=0 - 0.218 ± 0.002. (15) 



3 Virtual and soft real photon emission 

Taking into account the accuracy level of 0.1% for determination of p/ po we will drop 

terms of order r e . We will distinguish 3 kinds of contributions to 5: from emission of 

virtual, soft real, and hard real photons in the rest frame of kaon: 5 = 8y + $s + 8h ■ 

Standard calculation (see Appendix A for details) allows one to obtain the following 

contributions: 

from the soft real photons 



2Ae 1 
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*s = -|(Le-2)ln— + -L e --Li + l-—^(l + 0(r e )) , (16) 

from the virtual photons dy = 8c + Splm that make up charged fermion's renormal- 
ization, 8c (throughout this paper we use Feynman gauge): 



8c 



2tt 



1 r 3 , 

■-L A + -lnr e 



In 



A 2 " 



La + In 



A 2 



3 



(17) 
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here La = ln(A 2 /M 2 ), A is ultraviolet momentum cutoff, the first term in the curly 
braces conies from positron, the second one from kaon; and for the diagram in figl(f) 
in the point like meson (PLM) approximation, 5plm- 



jplm : 
a 

~ 2^ 



1 o M n 

-L A - - In r e - 2L e + In — L e - 1 + 2 In 2 y + 2 Iny + 2Li 2 {\ - y) 

2 X A 



(18) 



When these contributions are grouped all together the dependence on A (fictitious 
"photon mass") disappears. According to Sirlin's prescription Q we set A = Mw- 
The result can be written in the form: 



l + Ss + Sc + Sp L M = S W [1 + 



1) In A 



TT 

~6~ 



Li 2 (i - y) - 77 ln y 



In A 

Sw = 1 



3a 

47T ' 



(19) 



In the above equations L e = 21ny + ln(l/r e ), and Lw — h^M^ /M 2 ); M\y is the 
mass of W^, A = Ae/£v e , and Ae is the maximal energy (in the rest frame of kaon) 
of a real soft photon. We imply Ae < M/2. For the details of eqs (|l|), (0), ©, 
and ( pj| ) see Appendix A. 

Contribution from soft photon emission from structure-dependent part (such as 
for example, interaction with resonances and intermediate W^ 1 ) is small, of the order 



a Ae 

7T M 



< 1 



and thus is also neglected. 



4 Hard photon emission. Structure function ap- 
proach 

Next we need to calculate contributions from hard photons. We have to distinguish 
between inner bremsstrahlung (IB) and the structure-dependent (SD) contributions: 
Sh = Sib + Si nt + Ssd, where Si„t is the interference term between the two. The 
terms Si n t and Ss d are considered in the framework of the chiral perturbation theory 
(ChPT) to the orders of (p 2 ) and (p 4 ) and find their contribution to be at the level 
of 0.2% (see appendix C). 

Sh = 2 ^ {y z) { {L e - 1) U(y, z) - a {y, z) (2 In A + | 

2a {y,z)\n h ^ K ^+5 h ^ d , (20) 



(i 



with S hard given below. 

Extracting the short-distances contributions in form of replacement C — > CSw it 
is useful to split 5 (see eq.(2))in the form 

6(y,z) = S L + S NL , (21) 

where Sl is the leading order contribution, it contains the 'large logarithm' L e and 
Snl is the non-leading contribution, it contains the rest of the terms. 

Sl contains terms from 5c, 5s, Splm, and the contribution from the collinear 
configuration of hard IB emission (in the collinear configuration the angle between 
the positron and the emitted photon is small). Sl turns out to be 

6l = P^*(V,z). (22) 
2TTa (y,z) 

First we note that the kinematics of hard photon emission does not coincide with the 
elastic process (Region D, the strictly allowed boundaries of the Dalitz plot). In hard 
photon emission an additional region in the y,z plane, namely y < b_{z) appears. 
The nature of this phenomenon is the same as the known phenomenon of the radiative 
tail in the process of hadron production at colliding e + e~ beams. 

The quantity ^{y, z) has a different form for Region D and outside it: 

*(y, z) = *>(y, z), z > c(y), 2^fF e < y < 1 - rv ; (23) 

and 

z) = 4-<(y, z),z < c{y),2^T e <y<l-^~. (24) 

z) = when y > 1 — ^/r^ . Functions ^<, ^> are studied in Appendix D. 
Snl contains contributions from Sc, Splm, from SD part of hard photons and 
from the interference term of SD and IB parts of hard radiation. 

a 

Snl = ~v{z,y) , (25) 

IT 



where 



%(y, Z ) = 6 hard + - 
IT IT 



3 7T 2 3 

4 - y - L *2(i - y) - g ln y - ln (00) - y)/v) 



(26) 



and for the case when the variables y, z are inside D region: 



chard 



Z 2 (y,z); (27) 



2Tra {y,z) 

b(z)-y 



Z 2 (y,z) = -2Rphot 1D (y,z) + Rphot 2 D(y,z) + J dxj(x,y,z). (28) 

o 

Explicit expressions for Rjphot\^, and J are given in appendix B. The Born value 
and the correction to Dalitz-plot distribution A(y,z) = S(y, z)ao(y, z) is illustrated 
in tables 1,2. 
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We see that the leading contribution from virtual and soft photon emission is 
associated with the so called (5-part of the evolution equation kernel: 

(X \leading a fT u [ a o(t, z) p (l) (V\dt 

(Sc + Ss + S PLM ) a = - {Le -l)J -^— )P& (29) 

where 

Pj 1} (<)= 6(1 -t)(^21nA+0 . (30) 

The contribution of hard photon kinematics in the leading order can be found with 
the method of quasireal electrons [To|| as a convolution of the Born approximation 
with the #-part of evolution equation kernel Pg(z): 



„ « (Le - 1) / ^Jf) f (31) 



where 



^( Z ) = ^(l- Z -A). (32) 

In such a way the whole leading order contribution can be expressed in terms of 
convolution of the width in the Born approximation with the whole kernel of the 
evolution equation: 



1 )(«) = |m(jj 1) («)+Pi l) ( 2 )) . (33) 

This approach can be extended to use nonsinglet structure functions D(t, y) ||: 

&(*) 

dT LO (y,z) = J j dT (t,z)D(^, L e ), t = x + y, (34) 

max[y,b— (z)] 

D(z, L) = 6(1 -z) + ^LPW(z) + 1 (g) 2 pW(*) + ... , 

i 

pW( z )= / ^p(i)(a;)p(i-i)(f) , i = 2,3,.... 
J x x 

z 

One can check the validity of the useful relation: 

l 

dzP (1 \z) = Q. (35) 



The above makes it easy to see that in the limit m e — > terms that contain m e do 
not contribute to the total width in correspondence with Kinoshita-Lee-Nauenberg 



S 



(KLN) theorem Jl2j as well as with results of E.Ginsberg 0]. Keeping in mind the 
representation 



b(z) 

*(»,*) = / J «o(M)P (1) (f ) 



(36) 



maa:[y J b_(z)] 

one can get convinced (see Appendix D) that the leading logarithmical contribution 
to the total width as well as one to the pion spectrum is zero due to: 



l+r„ fc(z) 



f dz J dyV(y,z) = 0. 



(37) 



Using the general properties of the evolution equations kernels, eq(35) one can see 
that KLN cancellation will take place in all orders of the perturbation theory. The 
spectra in Born approximation are (we omit terms 0{r e ) ~ 10~ 6 ): 
For pion 



ld£o 
C dz 



M*) > (38) 

2 b[z) 2 

1 + ^R(z)) J dya (y, z)=U+ ^(z)) ~ (z 2 - ^f 2 , 



b-(z) 



and for positron 
ft \ 1 dT ° 



fo(y) 



i + 1 (b±\ g£ ~ r *- y} + 1 (K\ 2 y 2 {i~r„ 

i\rj l-y 6\rJ (l-j 



yf 



fo(y) = y2il - r *- y)2 . (39) 

The corrected pion spectrum in the inclusive set-up of experiment when integrating 
over the whole region for y (0 < y < b(z)) have a form 4>o(z) + (a/ir)(f>i(z) with 



<f>i{z) = ( 1 + — R(z 



dy[^< {y, z)lny- ao(y, z) In 



Hz) - y 
b-(z) - y 



1 ~ 



b(z) 



z 2(y,z)]+ J dy[^ > {y,z)lny + a (y 1 z)Z 1 {y,z) + -Z 2 (y,z)] 
b-(z) 



(40) 



the quantities Zi,Z% explained in Appendix E. This function do not depend on 
ln(l/r e ). Pion spectrum in the exclusive set-up (y,z in the region D) will depend 
on L e . It's expression is given in the Appendix E. 



9 



Numerical estimation of pion spectrum is illustrated in figure (3,5). 
The inclusive positron spectrum with the correction of the lowest order is f(y) 
(a/ir)fi(y) with f(y) given above and: 



i . r , (. a. 



h{y) = -(i e -l)/(y)- J a (y,z)\l + ^R(z)\ ln((6(z) - y)/y)dz+ 

c (y) 

2 1+r,r 2 

I~ \-\^V- L ^{l-y))}{y) + \ J Z 2 (y,z)(l + ^-R(z)) dz+ 

c (v) 

6(1- ^-y) f dz(l + ^± R{z)f[{l/2)Z 2 ~ a (y,z) In M ~ V }, (41) 



with 



/(//> =>(//> +(^Wtf)+(^ ±> ) J2(y), (42) 



jo(y)= J j J dza (t,z)P^(^) = 
y c(t) 

(2ta l^^ + | )A(! , ) + r|l±^ ta ^i + 

y 2 2(1 - y) rv 

-L(l - t> - y)[l - 5rv - 2r^ + y(4 - 13rv) - 17y 2 } ; (43) 

explicit expressions for j\(y) and j 2 {y) are given in Appendix D. 

Numerical estimation of positron spectrum is illustrated in figure (4,6). 

One can check the fulfillment of KLN cancellation of singular in the limit m e — » 
terms for the total width correction: f r " I(y)dy = 0. The expression for jo(y) agree 
with A(2) from the 1966 year paper of Eduard Ginsberg 

We put here the general expression for differential width of hard photon emission 
which might be useful for construction of Monte Carlo simulation of real photon 
emission in K e ^\ 

a dx dO 

dl™'" = eU'o 

with 



dT har d _1 rp (44) 

' Z7r x ZTTao(y,z) 



2uo 2Ae Ae Ae 

x = — > = y — , — < 1 ; (45) 

M M E e E e y ' 

and d0 1 is an element of photon's solid angle. The quantity T is explained in Ap- 
pendix B. 

For soft photon emission we have 



' Z7T X ZTT 



4 Te . y 



(1 - (ieCe) 2 1 - [3 e C e 



x < y-^ ■ (46) 
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Integrating over angles within phase volume of hard photon we obtain the spectral 
distribution of radiative kaon decay: 



dT a 1 



dTodx 2Trao(y,z) 



Zt + + y y r (fr a+ (*+y) a )( £ «- 1 ) +a!a ) 

2 , N R(z) 

-a {y, z) - 2( — ■ y) + J(x, y, z) 

x x + y 



yA < x < b(z) - y, A=— ^<1. 

5 Discussion 

Structure-dependent contribution to emission of virtual photons (see Fig 1 d), e)) 
can be interpreted as a correction to the strong formfactor of Kit transition, /+(£). 
We assume that this formfactor can be extracted from experiment and thus do not 
consider it. The problem of calculation of RC to K e ^ and especially the formfactors 
in the framework of CHpT with virtual photons was considered in a recent paper Jl4[ . 

As in paper Q we assume a phenomenological form for the hadronic contribution 
to the K — 7r vertex, but here we use explicitly the dependence of the form factor in 
the form 

U(t) = f+(0)(l + ^R(z)j . (47) 

We assume that the effect of higher order ChPT as well as RC to the formfactors can 
be taken as a multiplicative scaling factor for / + (0), which we take from of a recent 
paper Q. 

We assume such an experiment in which only one positron in the final state is 
present, but place no limits on the number of photons. The ratio of the LO contribu- 
tions in the first order to the Born contribution is a few percent, and for the second 
order it is about 



< 0.03% . (48) 



aL e s 
~2tt~ 

Due to non-definite sign structure of the leading logarithm contribution (see eq(|22|)) 
there are regions in the kinematically allowed area where \^(y, z)\ is close to zero. In 
these regions the non-leading contributions become dominant. 

The contribution of the channels with more than 1 charged lepton in the final 
state as well as the vacuum polarization effects in higher orders may be taken into 
account by introducing the singlet contribution to the structure functions. The effect 
will be at the level of 0.03% and we omit them within the precision of our calculation. 

The contribution of the 0(p 4 ) terms || turns out to be small. Indeed, one can 
see that they are of the order 0(aL r g , (p/A c ) 2 ) < O(10- 2 %), A c = 4nF n « l.2GeV 
(F-k = 93MeV is the pion life time constant), where p is the characteristic momentum 
of a final particle in the given reaction, p 2 < M 2 /16 ~ F 2 . So the terms of the orders 
Q{p A ) and 0{p 6 ) can be omitted within the accuracy of O (f x 10~ 2 ) < O (l0~ 4 ). 
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Our main results are given in formulae (2,21,22,26-28) for Dalitz plot distribution ; 
(38-41) for pion and positron spectra;(47) for hard photon emission;(53) for the value 
\V U8 \, in the tables and figures. The accuracy of these formulas is determined by the 
following: 

1. we don't account higher order terms in PT, the ones of the order of (aL e /ir) n , n > 
2 which is smaller than 0.03% 

2. structure-dependent real hard photon emission contribution to RC we estimate 
to be at the level of 0.0005. 

3. higher order CHPT contributions to the structure dependent part are of the 
order 0.05% f§ | 

All the percentages are taken with respect to the Born width. All together we believe 
the accuracy of the results to be at the level of 0.01. So the final result of our 
calculation may be written in the form 

^ = (1 + 6(1 ±0.01)) (49) 

the terms on the RHS are given in 

Here is the list of improvements comparing with the older calculations [|l], ^ : 

1. we eliminate the ultra-violet cutoff dependence by choosing A = Mw, 

2. we describe the dependence on the lepton mass logarithm L e in all orders of the 
perturbation theory and explain why the correction to the total width does not 
depend on m e , 

3. we treat the strong interaction effects by the means of CHPT in its lowest order 
0(P 2 ) and show that the next order contribution is small, 

4. we give an explicit formulas for the total differential cross section and explicit 
results for corrections to the Dalitz plot and particle spectra that might be used 
in experimental analysis. 

In the papers of E.Ginsberg the structure-dependent emission was not considered. 
T.Bccherrawy, on the other hand, did include it, and this will give rise to differences 
in the Dalitz plot. In addition, Ginsberg used the proton mass as the momentum 
cutoff parameter. 

We do not consider the evolution of coupling constant effects in the regions of 
virtual photon momentum modulo s qua re from the quantities of order M 2 up to 
M§, which can be taken into account |TTf (and for details see jlij) by replacing the 
quantity Sw by the Sew = 1 + («/tt) ^(M^/M?) = 1.0232. Taking this replacement 
into account our result for the correction to the total width is 

— = 1 + 5 = 1.02, (50) 
To 
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which results in 



\V us f+{0)\ =0.214 ±0.002. 



(51) 



So the correction to the total width is +2% while Ginsberg's result is —0.45% and 
Bccherrawy's is —2%. Neither Ginsberg nor Becherrawy used the factor Sew, and 
this factor (1.023) accounts for most of the difference between Ginsberg's and our 
result. Electromagnetic corrections become negative and have an order of 10~ 3 . The 
effect of the SD part, which E.Ginsberg did not consider, is small, of the order of 



We use the value of formfactor /+(0) = 0.9842 ±0.0084 calculated in the paper (mJ 
in the framework of ChPT, including virtual photonic loops and terms of order 0(p 6 ) 
of ChPT. To avoid double counting we use the mesonic contribution to /™ es (0) = 
1.0002 ± 0.0022 and the p 6 terms one /+(Q)| p 6 = -0.016 ± 0.0008). Our final result 



In estimating the uncertainty we take into account the uncertainties arising from 
structure-dependent emission ±0.005, theoretical errors of order ±0.0003, experimen- 
tal error ±0.0022 and the ChPT error in p 6 terms 0.0008. 

In tables 1,2 we give corrections to the distributions in the Dalitz-plot dT/(dydz) ~ 
a>o{y,z) + (a/ir)A(y,z). 

In figures (3-6) we illustrate the corrections to the pion and positron spectra. 
Here we see qualitative agreement for the positron spectrum and disagreement with 
the pion spectrum obtained by E.Ginsberg. 
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0.1% 



is: 



\V US \ = 0.21715 ± 0.0055. 



(52) 
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z/y 


0.07 


0.15 


0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


1.025 


3.83 


4.37 


3.71 


2.01 


-0.21 


-2.47 


-4.31 


-5.11 


-3.81 


0.975 




3.76 


3.49 


2.07 


0.05 


-2.07 


-3.83 


-4.61 


-3.35 


0.925 






3.26 


2.13 


0.32 


-1.67 


-3.35 


-4.11 


-2.88 


0.875 






3.04 


2.18 


0.58 


-1.26 


-2.86 


-3.60 


-2.39 


0.825 








2.25 


0.85 


-0.86 


-2.37 


-3.08 


-1.88 


0.775 










1.14 


-0.41 


-1.83 


-2.51 


-1.28 


0.725 










1.39 


-0.04 


-1.39 


-2.03 


-0.72 


0.675 












0.38 


-0.86 


-1.48 




0.625 














-0.35 


-0.89 




0.580 
















-0.23 





Table 1. Correction to Dalitz plot distribution A(y,z) 
(2))- 



ao(y, z)5{y, z) x 10 3 (see eq 



z/y 


0.07 


0.15 


0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


1.025 


0.0084 


0.069 


0.126 


0.163 


0.181 


0.178 


0.156 


0.114 


0.051 


0.975 




0.026 


0.089 


0.131 


0.153 


0.156 


0.139 


0.101 


0.0437 


0.925 






0.051 


0.099 


0.126 


0.134 


0.121 


0.088 


0.036 


0.875 






0.014 


0.066 


0.099 


0.111 


0.104 


0.076 


0.029 


0.825 








0.0337 


0.071 


0.089 


0.086 


0.064 


0.021 


0.775 










0.043 


0.066 


0.069 


0.051 


0.014 


0.725 










0.016 


0.044 


0.051 


0.039 


0.006 


0.675 












0.021 


0.034 


0.026 




0.625 














0.016 


0.014 




0.580 
















0.003 





Table 2. Dalitz plot distribution in Born approximation ao(y, z). 



7 Appendix A 



Here we explain how to calculate 5s, 5c, 5plm and how to group them into eq(^9[). 

Contribution from emission of a soft real photon can be written in a standard 
form in terms of the classical currents: 

(53) 



5s 



47ra 



(27T)* 



2lo 



P 
p-q 



Pe 



Pe ■ q 



where A is fictitious mass of photon. We use the following formulas: 
1 f d 3 q ( p 

,P ■ <?, 



2vr / 2w 



In 



2Ae 



1 



14 




15 




Figure 6: Correction to positron spectrum, fi(y) (see (41)). 
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1 f d 3 q ( p e V 



2tt J 2uo \Pe'qJ \ X J 2 

1 [d?q 2{p-p e ) _ T ^2Ae^ tt 2 1 r . 



2tt J 2w (p • «)(p e • g) - ie ln J ~ IT ~ I^ 6 ' (54) 

From them we obtain the eq(|l6|) . 

Consider now radiative corrections that arise from emission of virtual photons 
(excluding SD virtual photons). 

Feynman graphs containing self-energy insertion to positron and kaon Green func- 
tions (Fig.l,b,c) can be taken into account by introducing the wave function renor- 
malization constants Z e and Zk'- Mq —> Mo^xZe) 1 ' 2 . We use the expression for 
Z e given in the textbooks |l6| ; the expression for Zk is given in the paper fl5fl . The 
result is cq(|l7|) ■ 

Now consider the Feynman graph in which a virtual photon is emitted by kaon 
and absorbed by positron or by W-boson in the intermediate state (Fig.l,d,e,f). This 
long distance contribution is calculated using a phcnomcnological model with point- 
like mesons as a relevant degrees of freedom. To calculate the contribution from the 
region \k\ 2 < A 2 (A is the ultra violet cutoff) we use the following expressions for loop 
momenta scalar, vector, and tensor integrals: 

Re [** h fc2 =/ //* j f55 ) 

J iir 2 (k 2 - \ 2 )((k - p) 2 - M 2 )((k - p e ) 2 - m 2 e ) ' ' 1 ' 

A standard calculation yields: 

-1 fl M 2 1 1 

/= ^{2 ln ^ 2 " Le + ln2y + ^ 2(1 " 2/) "4 ln2re j ; (56) 



yM 2 \ 1 - y * i e \ 1 - y 

V hi y 

J = L A + f— y - + 1 . 

where L\ = ln(A 2 /M 2 ) and we omitted terms of the order of 0(m 2 /M 2 ). As a result 
we obtain 

rf 4 fc (l/4)Spp v (p + p')(-p e + k)(2p - k) Pe (p + p') 4 
ztt 2 (k 2 - \ 2 ){{k - p) 2 - M 2 )((k - p e ) 2 - m 2 e ) a ol2/,2jx 

i In 2 r e - 2L e + ln 

Z A 

In a series of papers JtJ A.Sirlin has conducted a detailed analysis of UV behaviour 
of amplitudes of processes with hadrons in 1-loop level. He showed that they are UV 
finite (if considered on the quark level) , but the effective cutoff scale on loop momenta 
is of the order My/, For this reason we choose 

L A = ln M - 



L A -^ln 2 r e -2£ e + ln^L e -l + 21n 2 y + 21ny + 2L; 2 (l-y)}> . (57) 



M 2 

The sum 5s + &c + &plm yields cq(p9h 
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8 Appendix B 

The matrix element of the radiative K e 3 decay 

K+(p) -> TrV) + e+ (Pe) + V{P») + 7(9) 

with terms up to 0{p 2 ) in CHPT || g § @ has the form 

G 



M 



hard 



(58) 



(59) 



where 



Q 



hard 



Q 



IB 



(P + p') 



Ql + Ql D = Qi B + Ql D 

~ (-Pe ~q + m e )j^ p^ 
2p e q pq 
n SD _ p 



(60) 



where tensor R pv describes (see eq(4.17) in Q) the structure-dependent emission (fig 
2(c)). 

Rau = gnu - -^-^ • (61) 

pq 

Singular at x = %Peq —> terms which provide contribution containing large logarithm 
L e arises only from Q e ^ . To extract the corresponding terms we introduce four-vector 
v = (x/y)p e — q, and x — is the energy fraction of the photon (9). Note that v — > 
when x ~ * 0. Separating leading and non-leading terms and letting f+(t) = 1, i.e. 
neglecting form factor's momentum dependence, we obtain: 



Sh 



r/r 



hard 



dr 



27rao(y, z) J x 



dx f dO^ 



2tt 



T ,x>yA 



(62) 



where 



T =ir E l s w(Q?r 



(1 + lb)v{pe 



ya (x + y,z) 
x + y 



r y 2 + (x + y) 2 (l-p e )(x + y) 



y 2 (l - /3 e C e ) y(l - /? e C e ) 5 



j/a (x + y, z) 



+ P. (63) 



The quantity V contains some non-leading contributions from the IB part and the 
ones that arise from the structure-dependent part: 



1 A/ 2 ^M 2 -T + y 



{l-x-y)\+ — 



M 2 x + y 



M 2 x 2 



with 



T v = 



4M 4 



T v + -T lv ) - (T RR + 2T R ) , (64) 
q \ a; 



{t> + p')pv{p + ft')v 
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Trr = R^xR^a^pSp PulXPela] 



Tr = Rf.x^j^Sp p v (p + p ) 



P M (Pe + 9)7, 



Pel\ ■ (65) 



_pq x 

To calculate these traces we use the following expressions for the scalar products of 
the 4-momenta (in units M): 

P 2 = 1, q 2 = 0, pi = 0, p 2 = r n , p 2 e = 0, pp e = -, 

PP = |, M=|, PP^ = 2 ( 2 ~ y~ z ~ ^ ' 
p'Pv = - {I- x-y-r-* + A e ) , p'q = -(x - A e - A v ) , 

P'Pe = ^ (2/ ^ i? ( Z ) + ^) > PvV = \ A »i Pel = 

PePv = ~ (R(z) ~ A e - A v ) , pv = 0, p e u = -iyl e , 

= --A e , p'v = - A u + A 

2 y 2 V V 

p I/ v = -—(xA e + (x + y)Ajj , 

A„ = A v —R(z) . 

x + y 

Three terms in the rhs of ( |63| ) have a completely different behavior. 

The first one corresponds to the kinematic region of collinear emission, when pho- 
ton is emitted along positron's momentum. The relevant phase volume has essentially 
3-particles form: 

{d ^ r u = (^i^l^ pu§(p 2 )6 4 {p _ pe _ pv _ p/ _ CM = 

M 4 ^-f3 Tr zdzydyxdxd0 7 dC eTr S{l - x- y- z + r 7T + ^— t-^^(l - ^C en ) + = 
64 2/2 M z 

2 

V -M A —dO y xdxdydz . (66) 



x + y 32 

The limits of photon's energy fraction variation are yA < x < b(z) — y. The upper 
limit is imposed by the Born structure of the width in this kinematical region. 

The second term corresponds to the contribution from emission by kaon. The 
relevant kinematics is isotropic. 
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The kinematics of radiative kaon decay and the comparison of our and E.Ginsbergs 
approaches is given in Appendix F. 

The third term corresponds to the rest of the contributions which contain neither 
collinear nor infrared singularities. 

Performing the integration over photon's phase volume provided y, z are in the D 
region we obtain: 

Hz)-y 

dx f d0 7 f dx y 2 y 2 + (y + x) 2 x 2 
~T = / — - — ■ — —a (y + x,z)[ (L e — 1) H 



2ir J x (y + x) 2 V 2 V 2 

y A 

b(z)-y b( z )-V 

dx , 



X 
yA 



[a (y,z)+x(^--y)]+ [ dxj, (67) 
x + y J 



we obtain (f27n with 



&(*)-!/ 



Rphot 1D = f dx{^--y)=R{z)\J^--y{b{z)-y); (68) 
J x + y y 

o 



Hz)-V 



Rphot 2D (y,z)= dx Xa °. {y + X y f =- (R( z ) + y(2-z))ln^~ 
J [y + xy y 



^(b(z)-y)(2^4 + 4-2z-b(z) + y) , (69) 



and 



2 v w \ K z 

J(x,y,z) = - [ ^V. (70) 
xl 2tt 



One can check that the sum of RC arising from hard, soft and virtual photons do not 
depend on the auxiliary parameter A. 

We note that the leading contribution from hard part of photon spectra can be 
reproduced using the method of quasi real electrons . [] 

Now we concentrate on the contribution of the third term in the RHS of d63|). 

To perform the integration over the phase volume of final states it is convenient 
to use the following parameterization (see Appendix F): 

d 3 p'd 3 p e d 3 p v d 3 q 4 , tt 2 4 dC e dC n 

^ 4 = o/ofo o 6 (P~P -Pe-Pv- q) = P-n-r^M dydzxdx ? =— , (71) 

2e lhj e le v lio lb JD 



1 the formula (10) in [|lo[| should read 

dW^ , t (k) 



2e'd?o 0b 



=P3 + fc 
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with 

D = 0i{\ -C 2 -Cl- Cl + 2CC n C e ) , & = ^/l-^f , (72) 
C = cos(p e ,p'), C e = cos(q,p e ), C w = cos(g,p') . 
The neutrino on-mass shell (NMS) condition provides the relation 



1 - PnC = — 

yz 



x + y + z-l-r v -^(l-t3 v C„)-^(l-C e ) . (7?,) 



For the aim of further integration of V over angular variables we put it in the 
form: 

V = xPi 4^ + xPi + P 3 A e + P A A V + P b A v A e , (74) 

A e = ^(l-C e ) , 

xz 

A v = x — A e — (1 — PttCk) . 



and 



Pi = | (1 - x - y) ; (75) 

P 2 = —7-^- + I (z(2x + 3y + 1) + 2x 2 + Axy + 3y 2 - 2x - 3y - 2) ; 
x + y 2 

P 3 = 1 - z - y - ix (x + y + z) ; 

Pa = -1 + z + y + \xy ; P 5 = -1 • 
Angular integration can be performed explicitly, we have 
f PvdCit y f p v C v dCit y(x + y -yt) 

J ww = 71 'J = z&aw mz) - {x + y) {2 - z)+ xyt)] ' 

(76) 

with 

A = (x + yf - 2xyt, t = 1 - C e . (77) 
Performing the integration over C v we have: 

1 fdC w p w 2y ( ( z R \ 1 \ 



wVD A 3 / 2 V V 2 « + 2/ / 2 



A 1 / 2 
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The following integrals are helpful in integrating the above expression. We define 

f 2 dtV' 

L 7^ 



f 2 dtt m 

C= / -=, to = 0,1, 2, 3; n= 1,3. (79) 



Then 



T o _ 4 T i_ S(x + y + a) 
1 "a' 1 ~ 3<r 2 ■ 

If = ^ (3a 2 + 3(z + y)p + 5(x + y) 2 ) . 

n = ^— v 4 = A> / 3 2 = ^(2(x + 2/ ) + a). 
pcr(a; + y) pa z Apa 6 

32 

^1 = ( ct2 + + v)p + 4^ + vf) ■ 



(80) 



where p = \x — y\ and a = x + y + p. 

The first term in dT hard together with the leading contributions from virtual and 
soft real photons was given in the form required by RG approach eq(]3(]). 

The non-leading contributions, 5 hard from hard photon emission, include SD emis- 
sion, IB of point-like mesons as well as the interference terms. It is free from infrared 
and mass singularities and given above (^) with 

2 2 

J(x, y, z) = P 1 R 1 + P 2 yl° + P^—l\ + |p 4 i? 4 + ^P 5 R 5 , (81) 

and 

Ri = "I- (V - x) ((2 -z){x + y)- 2R(z)) J° - y 2 ((x + y)I° - xll) , 
x + y 

R 4 = (2- z) J° - yl{ + (2R(z) -(x + y)(2 - z))((x + y)/ 3 - yll) + 

xy{{x + y)ll - yll) , 

R 5 = (2-z) I\ - ylf + (2R(z) - (x + y)(2 - z))((x + y)l\ - yl 2 ) + 

xy({x + y)I 2 - yll) . 

9 Appendix C 

The contribution to 5 hard from SD emission have the form: 

b{z) 

S h S D d = 2lra ^ z) J dxJ SD (x, y, z) , (82) 
o 

where 

2 2 

J SD {x,y, z) = Q X R X + yQ 2 I° 1 + + + ^Q 5 R 5 , (83) 
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with Ri given in Appendix B and 

Qi = -\y{x + y) , 

Q 2 = \ [2x(x + 2y + z-2) + 3y(y + z - 2)] , 
Q 3 = ~ [-8 + (z + y){A + 3a;) - 2.x + 3a; 2 ] , 

Qi = \ [4y + 4a; + 3xy] , Q 5 = -| • (84) 
The contribution to the total width have a form: 

/ / dydz{\ + ^R(z)) 2 JdxJ SD (x, y, z) 

5 SD = °— x . (85) 

2tt ffdydza (y,z)(l + ±R(z)) 2 

Numerical estimation gives: 



S SD = -0.005. (86) 



10 Appendix D 

The function <3/, defined as 

b(z) 

*(y,z)= J ja Q {t,z)P^){ V - t ), (87) 

b-(z) 

contains a restriction on the domain of integration, namely t exceed y or equal to 
it,which is implied by the kernel P^\y/t). Explicit calculations give: 

b(z) 2 2 

*<(*,*) = / fao(M)|^ = TO- y (2-,)]ln^ + 

b-(z) 

2a (y, z) In ^~* y + \{b{zf b. (z) 2 ) , 

2V^<y<i-\Av, = 0,1/ > l - V^v • (88) 

and 

6(z) 



*>(l/,*) = | ya (t,^ (1) (f)=«o(y^)[21n 



K*)-y | 3, 
y 2 J 



i(fe(^) 2 - y 2 ) + (b(z) -y)(2-y-z + b.{z)) + [R(z) y(2 z)] In (89) 

2 y 
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One can convince the validity of the relations: 

c iv) l+r* 



jo(y)= J dz* < (y,z) + J dz^faz); (90) 

and 

b-(z) b(z) 

J dyf < (y,z)+ J dy® > (y,z) = 0. (91) 

6_( z ) 

The last relation demonstrates the KLN cancellation for the pion spectrum obtained 
by integration of the corrections over y in the interval < y < b(z). 
The explicit expressions for ji(y) and j2{y) are: (for jo(y) see (fl2|)). 

■ / \ 2/ 3 (l - ?V - y) 3 /„, 1- TV -v3 

[3(1 - y)(l + y 2 ) + r.(y 3 + 3y 2)] In l^L- 

l mi l n _ y )2 [C 1 - ( 4 V " " 3y - 1 + rv(832/ 2 + 26y + 11) + 3r 3 ) + 

r 2 (31y 3 -15y 2 -39y + 47)] , (92) 



12(1 



12(1 - yf 



In 



i-y 



r 2 (y 4 + 6y 2 -8y + 3) 



y 

6(1 + y 2 )(l - yf - 4tv(1 - y)(2y 3 - y 2 + 4y - 3) + 
1 ~ ^ ~ V -(1 - y) 3 (247y 4 - 88y 3 - 28y 2 - 8y - 3)- 



720(1 - y) 3 
rv(l - y) 3 (733y 3 + 341y 2 + 129y + 57)- 
r 2 (1 - y)(707y 4 - 808y 3 + 212y 2 - 408y + 717) 

+ r 3 (173y 4 - 72y 3 - 492y 2 + 1048y - 477) - 12r 4 (l - y) 3 



(93) 
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Collection of the relevant formulae. 

The Dalitz-plot distribution in the region D: 



dT 



CSew dydz 



1 + A, 



t 



a 1 1 

ao(y, z)-\ — [-(L e -l)*>(y, z)+a (y, z)Z 1 + -Z 2 

7T Z Z 



Zi = \ - y - \ my - ln((b(z) y)/y) Li 2 (l - y) . (94) 
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The function Z 2 is defined in (f28|). Correction to the total width (we include the 
contribution of the region outside the region D), T = Tq(1 + 5): 



l-TV 



1 + 6 = S 



EW i 



n J J dzdya (y, z) (l + j^R(zf) 



I(y) hiydy+ 



b-{z) 



dz(l + ^R{z)) 2 [ I dy[-a (y,z)hx 



Kg) - v 

b-(z) - y 

b(z) 



(l/2)Z 2 (y,z)}- 



dy[a (y,z)Z 1 + (l/2)Z 2 



b-(z) 



with 



Z 2 (y, z) = Rphot 2 A(y, z) - 2Rphot 1A (y, z) + 

b(z) 



RphotiA(y, z) = R{z) In 
Rphot 2A (y,z) = 



h{z) - y dxx 



b(z)-y 

dxj(x, y, z); 

b-(z)-y 

y(b(z)-b.(z)) 



b-(z)-y 



(x + y) 



: a a (x + y, z) 



(b(z) - b.(z))(l -*- + 2y)- (y(2 -z) + R(z)) In ^ 



(95) 



(96) 



The expression in big square brackets in right-hand side of (95) can be put in the 
form: 



4>\(z)dz 



h{y)dy 



-0.035 (97) 



which results in 8 — 0.02. For the aim of comparison with E. Ginsberg result we must 
put here 

A+=0, I{y)=jo(v), M W =M P , (98) 

as was mentioned above we have reasonable agreement with E. Ginsberg results. For 
the inclusive set-up of experiment (energy fraction of positron is not measured) we 
have for pion energy spectrum given above ©. When we restrict ourselves only by 
the region D the spectrum becomes dependent on ln(l/r e ): 

ck^Tz = {MZ)+ J( 1 /2)P(^)(ln(l/r e ) - 1) + 



b(z) 



1 



dy[^>{y, z)lny + a (y, z)Z\ + -Z 2 \ 



l + ^R(z)\ , (99) 



&_(*) 
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with 



P(z) = \b-{zf{W{z) + &_(*)) In £± + \{b{z) 6_ (z)) 3 In ^^M - 

1 



6 &_(*)(6(z) - 6_(z))(3Mz) + *(*))• (100) 
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Our approach to study the radiative kaon decay has an advantage compared to the 
one used by E. Ginsberg - it has a simple interpretation of electron mass singularities 
based on Drell-Yan picture. The || approach to study noncollinear kinematics is 
more transparent than ours one. We remind the reader of some some topics of [|| 
paper. One can introduce the missing mass square variable 

l=(p v + k) 2 /M 2 = A v = (M -K- E e ) 2 /M 2 - {p^+p e f/M 2 . (101) 

the limits of this quantity variation at fixed y, z are pu by the last term: for collinear 
or anticollinear kinematics of pion and positron 3-momenta. Being expressed in terms 
of y, z they are (we consider the general point of Dalitz-plot and omit positron mass 
dependence) : 

< I < b-(z)(b(z)-y), (102) 
for the y, z in the D region and 

b(z)(b-(z) -y)<l< b-(z)(b(z) - y), (103) 

for the case when they are in the region A outside D: 

< y < b^(z),2^ < z < 1 + rv. (104) 

For our approach with separating the case of soft and hard photon emission we must 
modify the lower bound for I in the region D. It can be done using another represen- 
tation of I: 

I = x[l - (y/2)(l - C e ) - (z/2)(l - pC«)], (105) 

with C' e , CV-the cosine of the angles between photon 3-momentum and positron and 
pion ones, [3 — y/1 — m 2 / E%-is the pion velocity. Maximum of this quantity is b(z). 
Taking this into account we obtain for the region oh hard photon 

x > 2Ae/M = yA, A = Ae/E e « 1, (106) 

for the region D: 

yA < x < b(z) - y, yb(z)A < I < b_{z)(b(z) - y); (107) 
and for region A: 

b-(z) - y < x < b(z)-y, b{z){b^{z) — y) < I < b-(z)(b(z) - y). (108) 
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In particular for the collinear case we must choice C e = 1; C w = — 1, which corresponds 
to x + y < b(z). Let infer this condition using the NMS condition : 

(P k - Pe - Pv - kfjM 2 = R(z) -x-v+ {xy/2)(l - C e )+ 

(xz/2)(l - 0C n ) + (yz/2)(l - /3C ev ) = 0. (109) 

In collinear case we have C e = 1; C n = C ew .From NMS condition we obtain 1 — /3C W = 
(2/z(x + y))(x + y — R{z)). Using this value we obtain l co u = R(z)x/(x + y). Using 
further the relation R(z) = b(z)b-(z) we obtain again x < b(z) — y in the case of 
emission along positron. 

Comparing the phase volumes in general case calculated in our approach with 
using NMS condition with j| approach we obtain the relation: 

d0 i _ f ji f j.. f^._ f d 3 p l ,d 3 kS i (P-p l/ -k) 



The non-leading contribution arising from hard photon emission considered above: 

hB = j~j^-V IB , (111) 

with 

Vib = xG^ + xG 2 + G 3 A e + G A A V + G 5 A e A„, 
A e 

Gi = f(2-<,-z), 

g, = ^ + £ + + 2 V) + l(2 Z + My + *)) - 1 ; 

G 3 = -^x 2 ~ ^x{2 + z + y)-~(y + z); G 4 = gx(4 + y) + ^y-l; 

G 5 = -~ (H2) 

(note that Gi + Qi — Pi, see appendix B and C) can be transformed to the form: 

(R(z)-l) 2 



7/ B = (l/4) / d/[4-2y-4z-(l/4) J R(z) + (l/4)Z + yln- 



I 

2 ln u,i R <o )2 ^ + + ^ 2 ^ y + z ) - 2 + (V4)K4 + v)]ho- 

l(l + y(2-z)) 

(l/2)/i_! - ((1/2)/ + y + 2)7 2 _a + JJ. (113) 
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Here we use the list of integrals obtained in the paper of || : 

1 



= / d~f 



{kP K /M 2 ) m {kp e /M 2 ) n 



I 



10 



s 2 



y — z + s 



y — z — s 
/-i.o = (2-y-«)/2; 

hi 



, = 4//; % = 1 
hl = yl ln T 

(R(z)-ir 



in 



h-i 



R{z) - I ' lr e 
R(z)(2-y-z)-(2 + y-z)l 



2l(y(2-y- z)-2R(z) + 2l) 2 
~ In ■ 



y — z + s 



2 — y — z — s 



h-i = 



2{y(2-y-z)+2l-2R(z)) 



R(z)(2-y-z)-(2 + y-z)l, 2-y-z + s 

3 ' 

s 2 — y — z — s 



Besides we need two additional ones: 

1 



= ^{2-y-zf-Al. 

° y 2 R{zf 



■In- 



(h, 



(k Pe /M 2 )(2{kP K /M 2 ) + y yR(z) 1(1 + y{2 - z))r e ' 

1 = 2l + ys + y(2-y- z) 

(kP K /M 2 )(2(kP K /M 2 ) + y ys U 21 + ys - y(2 - y - z) ' 



(114) 



(115) 



One can see the cancellation of mass singularities (terms containing ln(l/r e )) in the 
expression Ijb- 

Numerical calculations in agreement (within few percent) of this and the given 
above expressions. 
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